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We introduce a thermal dynamics for the Ising ferromagnet where the energy variations occurring 
within the system exhibit a diffusive character typical of thermalizing agents such as e.g. local- 
ized excitations. Time evolution is provided by a walker hopping across the sites of the underlying 
lattice according to local probabilities depending on the usual Boltzmann weight at a given tem- 
perature. Despite the canonical hopping probabilities the walker drives the system to a stationary 
state which is not reducible to the canonical equilibrium state in a trivial way. The system still 
exhibits a magnetic phase transition occurring at a finite value of the temperature larger than the 
canonical one. The dependence of the model on the density of walkers realizing the dynamics is also 
discussed. Interestingly the differences between the stationary state and the Boltzmann equilibrium 
state decrease with increasing number of walkers. 



PACS Numbers: 2.70.Uu, 2.70.Tt, 5.50.-hq, S.lO.Ln 

I. INTRODUCTION AND MOTIVATIONS 

Consider an Ising ferromagnet consisting of an assem- 
bly of N spins, each placed at a site of a d-dimensional 
Euclidean lattice. Let lowercase italic letters, i, j,- ■ ■, de- 
note such sites, and Si, Sj, ■ ■ the relevant spin variables, 
so that the energy pertaining to a given configuration 
s = {sili^i has the form 

1 ^ 

i—l jr^i 

where the symbol ~ in the second sum restricts it to the 
sites j adjacent to i. P 

As it is well known, the Ising modelEl is amenableL-of 
an exact solution on simple two dimensional latticesoB, 
whereas for more complex lattices or higher dimensions 
the study of its thermodynamic properties strongly relies 
on numeric simulations. According to the so-called dy- 
namic Monte Carlo method, the thermodynamic canon- 
ical average of a generic observable ^'^(s) is obtained as 
a simple algebraic average over the Markov chain of cop.- 
figurations produced by a suitable stochastic algorithmQ. 
Indeed each configuration is obtained from the previ- 
ous one so that, in the asymptotic limit, the probabil- 
ity Pt(s) that a given configuration s occurs at tem- 
perature T is proportional to the Boltzmann canonical 
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factor, exp [— T^^£'(s)], independent of the initial con- 
figuration. 

One of the most commonly exploited algorithms is 
the Glauber's single-spin-flip algorithmtl. According to 
Glauber's prescription two subsequent configurations can 
differ at most by the value of a single spin variable. More 
precisely the spin of the configuration s relevant to the 
site k can flip with probability 

^^^^'^)^ l + cxp[T'^Ai...(s),] 

where 

A^fc(s) = 2sfc^s, (3) 

is the energy variation consequent to the process. 
Glauber's dynamics is completely deflned by Eq. (||) to- 
gether with a prescription for updating the spin system. 
The latter is commonly chosen in view of a computa- 
tional optimization and typically consists in an unphysi- 
cal sweep along parallel lattice linesQ. 

As we mentioned above such a dynamics has a com- 
putational rather than physical origin, having been de- 
vised as an efficient way to give estimates of the canon- 
ical averages. JMevertheless it is usually given a physical 
interpretationu: the energy variation consequent to the 
spin-fiips realizing the evolution of the system occur due 
to the coupling of the spin degrees of freedom to a heat- 
bath at temperature T. 

In the last decade a numbfjij^nf new dynamics for the 
Ising model were introducedtfuiJ, where the energy vari- 
ations are typically required to occur uniformly through- 
out the whole sample. 

Here we introduce an alternative dynamics exhibiting 
the diffusive behaviour of a random walk in an evolving 
landscape generated by the energy of the spin couplings 
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in the real space. The spin- flips are induced by a ran- 
dom walker hopping across the sites of the underlying 
discrete structure. The motion of such a random walker 
is affected by the spin interaction, being biased towards 
those sites where a spin-flip is energetically more favor- 
able. It should be remarked that a relaxation dynamics 
occurring under ihe action of random walkers was intro- 
duced by CreutzEEl, and subsequently exploited by many 
authors. However our walkers act quite differently from 
Creutz's demons. Indeed the latter diffuse freely, unbi- 
ased by the magnetic configuration of the system, and, 
in order to simulate the micro-canonical ensemble, they 
induce a spin flip only if they can afford it, according to 
the individual energy they are endowed with. 

In fact our model in inspired by the physical behaviour 
of manganitesllil, such as La2:Cai_a:Mn03, where the spin 
dynamics is determined by the presence of diffusing exci- 
tations interacting with the magnetic degrees of freedom. 
In particular, for the manganites the excitations are given 
by the charged carriers, electrons and vacancies, present 
in the non-stoichiometrical compounds. Obviously, this 
model is not intended to be a phenomenological descrip- 
tion of such complex systems, but it aims to evidence the 
influence of diffusing excitation interacting with spins on 
the thermodynamic behaviour of the paradigmatic model 
of all magnetic phase transitions. 

The analysis of the diffusive dynamics is carried out 
from the numerical point of view. All the simulations 
illustrated in the present paper refer to two-dimensional 
arrays of spins, so that the observed features of the sta- 
tistical model can be compared to the analytically known 
results pertaining to the canonical Ising model. 

The plan of the paper is as follows: in Section || we de- 
scribe and briefly discuss the diffusive algorithm. Section 
III and IV are devoted to the analysis of the limit situ- 



ation where the dynamics is realized by a single walker. 
In the former we provide numerical evidences that the 
magnetic system is driven to a thermodynamically well- 
behaved steady state which differs from the canonical 
equilibrium state in a non trivial-manner. Since a mag- 
netic phase transition is still observed, the consequent 
critical behaviour is analyzed in Section IV. The es- 



timated critical exponents do not deviate significantly 
from the values pertaining to the canonical Ising model. 
In Section ^ the results of simulations where the system 
evolves subject to more than a single walker are analyzed. 
Section consists in our conclusions. 



Hence, for any configuration s of the Ising system and 
for any value of the temperature T, the diffusion of the 
walker is biased towards those sites where a spin-flip is 
more likely to occur according to Glauber's probability, 
Eq. (^. Note that Equation (^) implicitly yields a prob- 
ability 



PT(s,i,«) 



1 



{2d)-' Y.p^{s,k) 



(5) 



that the walker does not move. If this is the case the mag- 
netic configuration of the system s remains unchanged as 
well. According to equations and the evolution of 
the probability P(s, i, t) that, at time step t, the mag- 
netic system is in the configuration s while the walker is 
located at site i, is governed by the master equation 

PT{s,i,t+l) - PT{s,i,t) = 

(2d)-i^[FT(s:;,j,t)p?«,z)+ 



(6) 



where sj, is a shorthand notation for the configuration 
obtained flipping the spin Sk of s. Unfortunately the 
master equation approach is of very little use in ana- 
lyzing the long-time behaviour of this process. Actu- 
ally in the canonical heat-bath case the master equation 
was exploited to build a dynamics driving the system to 
a known asymptotic probability. This was obtained by 
imposing the quite restrictive detailed balance condition. 
Conversely, we introduced an evolution algorithm and 
our aim is the study of the resulting asymptotic state of 
the Ising system, if any. Hence we carry out our anal- 
ysis mainly by means of numerical simulations. For a 
better comparison we focus on the two-dimensional sys- 
tem, where most of the results pertaining to the canonical 
Ising model are exactly known. Hence all of our simula- 
tions refer to an Ising system consisting oi N = spins 
placed at the sites of a two-dimensional square lattice. In 
order to avoid time consuming procedures dealing with 
the walkers bumping into the borders of the system we 
adopt periodic boundary conditions. 

In the next two Sections we focus on the situation 
where the evolution of the Ising system is realized by 
a single walker. This is a limit case, in that, for a reason- 
ably macroscopic system, it corresponds to a vanishing 
density of walkers. The resultspertaining to larger den- 
sities are discussed in Section 0. 



II. DIFFUSIVE THERMAL DYNAMICS 

As we mentioned in the previous section the relaxation 
dynamics we propose is realized by means of random 
walkers diffusing through the sites of the Ising system. 
The probability that a walker located at site i hops on 
site J and flips the relevant spin, where j is one of the 2d 
neighbours of i, is given by 

PT{s,i,j)^{2d)-'p^{s,j) (4) 



III. NON-CANONICAL EQUILIBRIUM STATES 

In the following we give numerical evidence that the 
spin-flip dynamics realized by a single walker drives the 
magnetic system to a thermodynamically well-behaved 
steady state which differs from the canonical equilibrium 
state in a non trivial manner. All the simulations we 
performed show that, at a given value of the external pa- 
rameter T the system eventually reaches a steady state. 



3 



0.36 



0.35 



0.34 



0.33 



0.32 



0.31 

















1 



100 



200 



400 



800 



1600 



FIG. 1: Finite size scaling for the specific energy of an Ising 
system subject to the diffusive dynamics described in Section 
O at T = 2.4. All the measurements were carried out in the 
stationary regime. The error-bars represent the fluctuations 
about the average values. 
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FIG. 2: Finite size scaling for the specific magnetization of 
an Ising system subject to the diffusive dynamics described 
in Section ^ at T = 2.4. All the measurements were carried 
out in the stationary regime. The error-bars represent the 
fluctuations about the average values. Note that despite the 
symmetry of the Hamiltonian the measured value for the spe- 
cific magnetization is non-zero. Similar to the canonical case 
the Ising system displays a phase transition. 



characterized by a well defined value of the (time) aver- 
age of the macroscopic observables. These features are 
clearly recognizable in Figures |l| and ^, where the average 
values of the specific energy e(s) = N~^E{s) and mag- 
netization m(s) = Si=i plotted for systems 



with different sizes, at a fixed value of the temperature 
parameter. Furthermore, as one would require, the entire 
sample and any reasonably macroscopic portion of it are 
characterized by the same specific values of the thermo- 
dynamic observables. As it is clearly shown in Figures 
^ and ^, the fluctuations about the average value of the 
specific thermodynamic observables exhibit the expected 
scaling behaviour, decreasing as the inverse square root 
of the system size. Figure ^ shows the average values of 
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FIG. 3: Finite size scaling for the fluctuation about the aver- 
age value of the specific energy for an Isin g sy stem subject to 
the diffusive dynamics described in Section at T = 2.4. The 
slope of the linear fit (dotted line) of the measured data (•), 
—0.48 ± 0.02, is in good agreement with the expected value, 
0.5. All the measurements were carried out in the stationary 
regime. 

the specific energy and magnetization at different values 
of the temperature parameter, T = Note that, sim- 
ilar to the canonical Ising model, the system exhibits a 
magnetic phase transition. We verified that the long time 
behaviour of the macroscopic observables of the system 
is uniquely determined by the only external parameter 
characterizing the dynamics, namely the temperature T. 
In particular different choices for the initial configuration 
of the magnetic system do not affect the asymptotic be- 
haviour of the system, except for a possible trivial bias on 
the orientation of the spontaneous magnetization. Sim- 
ilar to what happens in the heat-bath case, both of the 
possible orientation of the specific magnetization are in 
principle equally likely, but a strongly magnetized initial 
configuration is very likely to evolve in a stationary state 
exhibiting a spontaneous magnetization along the same 
direction. When the initial configuration has no net mag- 
netization there is no bias on the direction of the sponta- 
neous magnetization exhibited by the system below the 
critical temperature. In this situations the formation of 
macroscopic domains exhibiting opposite net magneti- 
zation is observed (see Fig.^. In the long time regime 
one of the domains eventually prevails against the others. 
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FIG. 4: Finite size scaling for the fluctuation about the av- 
erage value of the specific magnetization for an Ising system 
subject to the diffusive dynamics described in Section |^ at 
T — 2.4. The slope of the linear fit (dotted line) of the mea- 
sured data (•), —0.48 ± 0.02, is in good agreement with the 
expected value, 0.5. All the measurements were carried out 
in the stationary regime. 
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FIG. 6: Typical configurations for a 400 x 400 Ising array 
with magnetization (m) = 0.86 for the diffusive dynamics (left 
panel) and the Glauber's dynamics (right panel). Note that 
the same magnetization is attained for different temperatures. 
More precisely T"^'^ = 2.4 > r'='"s and t'^^^'"'"''' = 2.13 < 
T'"'°8^ ^here TJ''"«5 = 2/ln(l + V2) ~ 2.269 is the critical 
temperature of the canonical Ising model. Note further that in 
the diffusive case the domains present smoother boundaries. 



system evolving under the action of a single walker. How- 
ever, such a quantitative difference cannot be accounted 
for by means of a simple rescahng of the temperature. 
Indeed if the diffusive dynamics acted as heat-bath dy- 
namics with a rescaled temperature, T' = t^T), the 
generic configuration s would occur with a probability 
P(s,T) cx exp — r' ^E{s) in the asymptotic regimeQ. 

Hence the joint probability for the occurrence of a con- 
figuration such that e(s) = e' and m(s) ~ m! would be 
of the form 
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V{e,m\T) = Z(T)J(e',m')e 



(7) 
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where 



J(e', m') = 5{e' - e(s)) 5{m! - m(s)) (8) 
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and 



Z{T) = Y,J{e',m') 



-T'^^Ne' 



e ^ - =^e-^'-^^(^) (9) 



FIG. 5: Macroscopic observables for a 400 x 400 Ising ar- 
ray. Filled circles (•): specific magnetization; open circles 



(o): specific energy. 



As we already mentioned the system exhibits a magnetic 
phase transition. From Fig. ^ and ^, it is clear that the 
critical temperature for the diffusive dynamics is signifi- 
cantly larger than the value pertaining to the canonical 
Ising model. Such value, T^^^s = 2/ln(l + V2) ^ 2.269 
is analytically known and the heat-bath dynamics yield 
quite satisfactory numerical estimates of it. As we will 
see in the next Section, accurate estimates yield the value 
Tc = 2.612 ±0.002 for the critical temperature of an Ising 



Last equality of Eq. ^ was obtained making use of Eq. 
(^). Now note that, according to Eqs. (Q) and (||), 



^(e,m,T) 



P(e,m,T) 



J(e, m) 



(10) 



Hence, if the rescaling hypothesis was true, the plot of 
V{e,m,T) versus m at a fixed value of the specific en- 
ergy would not depend on the temperature: the curves 
pertaining to the same specific energy at different val- 
ues of the temperature T would overlap. This can be 
verified TCit|h great precision in the case of the heat-bath 
dynamicsE3 The numeric estimates of these curves for two 
different values of the specific energy e are plotted in fig- 
ures and ||. The fact that the curves pertaining to the 
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FIG. 7: Normalized joint distributions V(e, m, T) for a 20 x 20 
Ising array subject to the diffusive dynamics. All of the curves 
refer to the specific energy e = 0.220. Filled circles (•): T = 
2.25; open circles (o): T = 2.50; crosses (x): T = 2.75; dotted 
line: heat-bath dynamics (temperature independent). 




FIG. 8: Normalized joint distributions V{e, m, T) for a 20 x 20 
Ising array subject to the diffusive dynamics. All of the curves 
refer to the specific energy e = 0.345. Filled circles (•): T = 
2.50; open circles (o): T = 2.75; crosses (x): T = 3.00; dotted 
line: heat-bath dynamics (temperature independent). 



diffusive dynamics at different temperatures are well dis- 
tinct from one another and from the (temperature inde- 
pendent) curve characterizing the canonical Ising model 
proves that the differences between the diffusive and the 
heat-bath dynamics cannot be accounted for by means of 
a simple rescaling of the temperature. 

In summary the steady state asymptotically reached 
by an Ising system subject to the action of the diffusive 
dynamics is thermodynamically well-behaved, and yet it 



is non-trivially different from the canonical equilibrium 
of the Ising model. Hence we will refer to the steady state 
of the diffusive dynamics as a non-canonical equilibrium 
state. 

In the following section we analyze the critical be- 
haviour of the Ising system subject to the diffusive dy- 
namics. 



IV. CRITICAL BEHAVIOUR 

As we already observed in the previous Section, an 
Ising system subject to a single diffusive walker exhibits 
a magnetic phase transition at a finite value of the tem- 
perature (see Figure ^). Similar to what happens with 
the Ising model, this phenomenon is accompanied by a 
singular behaviour of the thermodynamic functions. In 
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FIG. 9: Critical exponent, f3, of the magnetization for an Ising 
system subject to the diffusive dynamics. The measures were 
performed on a system consisting of A*' = 2000 x 2000 spins. 
The error bars denote the standard deviation of the measures. 
The dotted line is the best y = A jT-Td''. The estimated 
values for the critical temperature and for the exponent are 
Tc = 2.612 ± 0.001 and /3 = 0.127 ± 0.002, respectively. The 
latter is consistent with the relevant critical exponent of the 
canonical Ising model, /Rising = 1/8. The vertical dashed line 
indicates the critical temperature. 



Figure |^ accurate data for the magnetization of an Ising 
system subject to the diffusive dynamics are plotted. The 
error bars represent the standard deviations about the 
average values. Similar to the canonical case, the data 
show a good agreement with a critical behaviour of the 
form m(T) ~ \T — Tc\^ , where the estimated values are 
for the critical temperature and for the critical exponent 
are Tc = 2.612 ±0.001 and /? = 0.127±0.02, respectively. 
The latter result is in good agreement with the relevant 
critical exponent of the two-dimensional canonical Ising 
model, /?ising = 1/8. As we already discussed in the previ- 
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ous Section, the critical temperature is appreciably larger 
than the canonical value, T^^""s « 2.269. 

Figure ^ shows the specific heat c(T) = d(e)T/dT 
for a 400 x 400 Ising array subject to the diffusive dy- 
namics as a function of the temperature T. The vertical 
dashed lines are placed at the critical value of the temper- 
ature. The dotted curves refer to functions of the form 
/(T) ^ a + b log(|T - Tel), and they fit quite satisfacto- 
rily the data. Hence, similar to the canonical case, there 
is a strong signature of a logarithmic divergence of the 
specific heat at the critical temperature. Figure pi] shows 
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FIG. 10: Specific lieat for an Ising system subject to diffusive 
dynamics (•). The dotted curves fitting the data are of the 
form /(T) = a -I- 6 log(|r - Td). The vertical dashed line 
indicates the estimated value of the critical temperature, Tc ~ 
2.612 ± 0.002. The diamonds (o) refer to the quantity Cf (T) 
defined in Eq. @. 



a log-log scale plot of the quantity 

N 



XiiT) = - [{m')T - {m)l 



(11) 



versus \T — Tc\. The data are consistent with a power 
law of the form Xf(r) - |T - Td'^, and the slope 
7 = 1.73 ± 0.06 of the linear fit is in good agreement 
with the critical exponent governing the behaviour of the 
same quantity in the of the canonical case, 7ising = 7/4. 
We recall that for a system at canonical equilibrium 
the fluctuation-dissipation theorem yields the relations 
c(r) = Cf(T) and x{T) - Xf{T), where 



Cf{T) 



N 

y2 



(12) 



and x(.T) — d {m)T.h/d h is the magnetic susceptibility 
of the system. In our case this is not necessarily true, 
since the equilibrium distribution for a diffusive dynam- 
ics at the temperature T is not proportional to the canon- 
ical Boltzmann factor exp as we discussed in 



the previous section. Hence the critical behaviour of two 
quantities canonically related by a fluctuation-dissipation 
relation should be explicitly analyzed and compared. A 
preliminary step in this sense is shown in figure |lO|, where 
the diamonds represent the estimates of the quantity 
Cf(T) defined in Eq. (|l^). These results seem to indicate 
that the diffusive character of the dynamics does not pro- 
duce large deviations from the fluctuation-dissipation re- 
lation. Of course this result needs to be checked through 
the analysis of the response of the system to an external 
magnetic field. However we mention that other studiesE^ 
of Ising systems subject to dynamics which do not yield 
the canonical equilibrium state actually estimate the crit- 
ical exponent 7 in terms of the quantity XtC^) defined in 

Eq. dnl). 
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FIG. 11: Log-log scale plot of the magnetic susceptibility Yf 
(defined in terms of the magnetic fiuctuations, see Eq. ([ll|)) 
versus |T — Tel. The estimated slope of the straight line fit- 
ting the data is —1.73 ± 0.07. This value is consistent with 
the critical behaviour of the magnetic susceptibility for the 
two-dimensional canonical Ising model, characterized by the 
critical exponent 7ising = 7/4. 



For temperatures sufficiently far from the critical value 
the correlation function 



rT(i, j) = {Si Sj)t - {Si)t {Sj)t 



(13) 



exhibits an exponential decrease as the distance between 
the relevant sites , rij = \ \i — j\ 



TxiiJ) (X exp 



mcreases: 



1 1] 



(14) 



As it is shown in Fig. 12 the correlation distance ^ (T) fea- 
tures a maximum very close to the critical temperature. 
The fact that this maximum suggests a critical tempera- 
ture slightly larger than the value previously (and more 
accurately) estimated must be ascribed to finite size ef- 
fects. We further mention that the correlation lengths 
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are roughly consistent with a power-law divergence of the 
form £,{T) ~ |T — Tc\~'^, the estimated exponent being 
remarkably close to the relevant critical exponent of the 
canonical two-dimensional Ising model, falsing = 1- As it 
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FIG. 12: Correlation length as a function of the temperature 
for a 600 x 600 Ising array subject to the diffusive dynamics. 
Due to finite size effects, the critical behaviour is exhibited for 
temperature slightly higher than the estimated critical value, 
Tc = 2.612 ± 0.002. 

is well known0, for temperatures very close to the crit- 
ical value, and hence for correlation lengths very large 
compared to the lattice constant, the correlation func- 
tion of the canonical Ising model is characterized by a 
power-law behaviour of the form TT{i,j) rfj. The rel- 
evant critical exponent is defined to he j] — p — d + 2, 
where d is the dimensionality of the Ising array. For the 
two-dimensional canonical Ising model its value is exactly 
known to be ryismg = 1/4. Such behaviour is expected to 
be observable in the distance range 1 ^ <C ^{T). As 
it is shown in Figure |l^ the data yielding the largest cor- 
relation length are not consistent with a pure power law 
behaviour, not even for relatively small distances. This 
suggests that the relevant temperatures, T = 2.63 is not 
sufficiently close to the critical value. Moreover it should 
be recalled that the finite size of the system can intro- 
duce appreciable deviations from the theoretical predic- 
tion, which is strictly true in the thermodynamic limit. 
Nevertheless we mention that a function of the form 

/(r)=aexp(-^)r-^ (15) 

fits the data quite satisfactorily over a wide range of dis- 
tances. The estimate for the exponent, rj = 0.252 ±0.002 
is remarkably close to the critical exponent T^ising of the 
two-dimensional Ising model. A fit of the same kind 
satisfactorily applies to the data relevant to the tem- 
perature T = 2.64, and it gives a critical exponent 
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FIG. 13: Log-log scale plot of the correlation function 
rT=2.63(*,i) versus the inter-spin distance, rtj. The data 
appreciably deviate from a power law behaviour. The dotted 
line refer to a function of the form defined by Eq. (p^. The 
exponent rj is quite close to the critical exponent jyising = 1/4. 
More precisely we estimate 772.63 = 0.252 ± 0.002. The same 
fit satisfactorily applies to the data relevant to T = 2.64, 
where the exponent r;2.64 = 0.266 ± 0.003 is once again in 
good agreement with the canonical value. 



r] = 0.266 ± 0.002 once again in good agreement with 
the canonical value. These results could be an indica- 
tion of the fact that, in the suitable regime, the corre- 
lation function TxaiT^ii j) actually exhibits a power law 
behaviour characterized by the same critical exponent as 
the two-dimensional canonical Ising model. 

The results illustrated so far suggest that a two- 
dimensional Ising ferromagnet subject to the diffusive 
dynamics exhibits a critical behaviour belonging to the 
same universality class as the canonical case. 

Before concluding the present section we mention a 
further result which has no counterpart in the canonical 
case. Indeed, in the diffusive case, it is possible to relate 
the configuration of the magnetic system to the position 
of the random walker causing its evolution. More pre- 
cisely it is possible to measure to what extent the pres- 
ence of the random walker at a given site i and the local 
magnetization at a given site j infiuence each other. Nu- 
merical estimates of such correlation, which we illustrate 
elsewheretJ, show that it exhibits a critical behaviour in 
correspondence of the critical temperature. 



V. RECOVERING THE CANONICAL 
EQUILIBRIUM: DENSITY OF WALKERS 

The results illustrated in the previous two Sections re- 
fer to the case of an Ising system subject to the action of a 
single walker. Here we analyze the results of simulations 
where the evolution of a square Ising array consisting of 
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N spins is realized by n non interacting walkers. The 
simulation is initialized giving each walker a randomly 
chosen position. Two or more walkers are allowed to oc- 
cupy the same site of the system. Actually this is very 
likely to happen for a sufficiently large density of walkers 
p — n/N . An elementary step of the simulation consists 
in the application of the algorithm described in Section 
H to all of the walkers, according to a fixed progressive 
sequence. When its turn comes, a given walker acts as 
if the other ones were not present. A qualitative argu- 
ment leads us to expect that in the case of an infinite 
density of walkers the results pertaining to the canonical 
Ising model are recovered. Consider a very large density 
of walkers, such that n ^ N. In this situation the ef- 
fect of a single step of the simulation is very similar to 
the action of the Glauber's dynamics with a random se- 
quence update. Indeed the probability that a walker fiips 
a given spin Sj is given by the product of the probabil- 
ity A/N that it is located at one of the four neighbours 
of site i times the hopping probability, p'^{s,i)/4:. But 
p'^{s,i)/N is exactly the probability that the spin Si is 
flipped at any step of the Glauber's heat-bath dynamics 
with random site update. Hence, in the limit case of an 
infinite density of walkers, the diffusive dynamics should 
necessarily drive the system to the same asymptotic state 
as the heat-bath dynamics, thus reproducing the results 
of the canonical Ising model. The same should be clearly 
true for a sufficiently large density of walkers, provided 
that the consequent number of elementary steps of the 
heat-bath dynamics drive the system sufficiently close to 
the canonical equilibrium state. The numerical simula- 
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FIG. 14: Critical temperatures with increasing density of 
walkers for a 400 x 400 Ising array. The horizontal dotted 
lines are placed at the critical temperatures of the diffusive 
dynamics with a single walker (Tc ~ 2.612) and at the critical 
temperature of the canonical Ising model (T^^'''"^ ~ 2.269). 



tions we performed evidenced that the diffusive dynam- 
ics drives the Ising system to a thermodynamically well 
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behaved stationary state, so that the average values of 
the specific energy and magnetization are characterized 
only by the temperature, T, and the density of walkers, 
p. In particular they confirmed the qualitative argument 
discussed above: the results pertaining to the canoni- 
cal Ising model at a given temperature T were recovered 
for a sufficiently large density of walkers. More gener- 
ally we verified that the average value of a macroscopic 
observable relevant to a finite density lies within the in- 
terval between the values pertaining to the limit situa- 
tions, namely the vanishing density limit p = N^^ ~ 
analyzed in Sections |l| and IV , and the canonical Ising 
model. The influence of the density of walkers is syn- 
thesized in Figure where the estimated values of the 
critical temperature are plotted. Note that the larger 
variation occurs for densities between 10~^ and 1. 



VI. CONCLUSIONS AND PERSPECTIVES 

In this paper we exploited the basic microscopic model 
of diffusion, namely the random walk, to build a dynam- 
ics for the Ising ferromagnet which retains the diffusive 
character of the thermal motion of a localized excitation. 
In this framework the density of walkers is related both 
to the intensity and spatial extent of the coupling of the 
spins to the thermal degrees of freedom of the system, 
and hence to the velocity of the thermal dynamics. The 
larger is the density of walkers, the faster is the dynamics 
and the closer to the uniform heath-bath are its effects. 

Even if our model is rather simplified with respect 
to the real physical systems which inspired it, it ex- 
hibits new important qualitative features which should 
be relevant even in more complex and realistic situations. 
Therefore it would be of great interest to verify whether 
the diffusive thermal dynamics gives a satisfactory de- 
scription for the slow evolution of a magnetic system with 
hopping excitations. This would of course require the 
extension of the study of its effect on the more realistic 
three dimensional lattice. In view of recognizing the ac- 
tion of a diffusive dynamics in a real system, the results 
of Section IV suggest that the study of the critical be- 
haviour is not particularly significant. Indeed the latter 
does not allow to discriminate among the canonical and 
the diffusive dynamics. 

On the other hand the results illustrated in section ^ 
indicate that the velocity of the non-equilibrium dynam- 
ics (related to the degree of coupling of the spins to the 
thermal vibrations, i.e. to the density of excitations in- 
duced in the system by an external source) significantly 
affects the critical temperature of the system. Hence the 
same system subject to diffusive thermal dynamics with 
different velocity should display different critical temper- 
atures. Furthermore we remark that Fig. || suggests 
that the velocity of the dynamics influences the shape 
and growth of the macroscopic magnetic domains. 

In particular a very slow diffusive dynamics gives rise 
to domains with smoother boundaries. This feature sug- 
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gests that the final shape of the magnetic domains within 
a physical system possibly governed by a diffusive ther- 
mal dynamics could be controlled by a fine tuning of 
the density of excitations mediating the coupling to the 
thermal degrees of freedom. A detailed study of domain 
growth in the diffusive dynamics will be the subject of a 
forthcoming papeuij. 
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